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The relation between saddle points of the potential of a classical many-particle system and the 
analyticity properties of its thermodynamic functions is studied. For hnite systems, each saddle 
point is found to cause a nonanalyticity in the Boltzmann entropy, and the functional form of this 
nonanalytic term is derived. For large systems, the order of the nonanalytic term increases un- 
boundedly, leading to an increasing differentiability of the entropy. Analyzing the contribution of 
the saddle points to the density of states in the thermodynamic limit, our results provide an expla- 
nation of how, and under which circumstances, saddle points of the potential energy landscape may 
(or may not) be at the origin of a phase transition in the thermodynamic limit. As an application, 
the puzzling observations by Risau-Gusman et al. on topological signatures of the spherical model 
are elucidated. 
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Phase transitions, like the boiling and evaporating of 
water at a certain temperature and pressure, are com- 
mon phenomena both in everyday life and in almost any 
branch of physics. Loosely speaking, a phase transition 
brings about a sudden change of the macroscopic proper- 
ties of a many-particle system while smoothly varying a 
parameter (the temperature or the pressure in the above 
example). The mathematical description of phase tran- 
sitions is conventionally based on (grand)canonical ther- 
modynamic functions, relating their loss of analyticity to 
the occurrence of a phase transition. Such a nonanalytic 
behavior in a (grand) canonical thermodynamic function 
can occur only in the thermodynamic limit in which the 
number of degrees of freedom TV of the system goes to 
infinity 

Many researchers took it for granted that the same 
were true also for microcanonical thermodynamic func- 
tions. Recently, however, it was observed that the mi- 
crocanonical entropy, or Boltzmann entropy, of a finite 
system is not necessarily analytic, and nonanalytic en- 
tropy functions of finite systems have been reported for 
certain classical models [2, H 0] • 

In light of their conceptual similarity to the definition 
of a phase transition, it is tempting to regard finite- 
system nonanalyticities of the entropy as phase transi- 
tionlike phenomena. This point of view is advocated in 
0| , and the authors of that reference argue that such non- 
analyticities should also be measurable experimentally, at 
least in very small systems. However, such an interpre- 
tation is complicated by the fact that, as discussed in 
for typical models the number of nonanalytic points of 
the entropy increases unboundedly with the number of 
degrees of freedom N. 



Because of their typically large number, one might as- 
sume that nonanalyticities of the finite-system entropy 
were unrelated to the occurrence of a phase transition in 
the thermodynamic limit. A theorem by Franzosi and 
Pettini [5| , however, indicates the existence of a relation 
between finite-system and infinite-system nonanalytici- 
ties (and we will come back to that theorem later). The 
purpose of the present Letter is to clarify and quantify 
the relation between saddle points of the potential en- 
ergy and nonanalyticities of the entropy of finite systems 
and, in a second step, to use this relation in order to ex- 
plain how, and under which conditions, a large number 
of "weak" nonanalyticities of the finite-system entropy 
may give rise to a single nonanalyticity (or a few) in 
the thermodynamic limit. At least in the common situ- 
ation of ensemble equivalence, such a nonanalyticity of 
the infinite-system entropy will correspond to a nonana- 
lyticity in the canonical free energy, and hence to a phase 
transition. 

Our approach will be the following: We start from the 
case of classical many-particle systems, announcing an 
exact expression for the leading nonanalytic term that 
each saddle point of the potential energy landscape con- 
tributes to the density of states. This result, apart from 
its essential role for further steps of our analysis, is of in- 
terest in its own, giving a model-independent and quanti- 
tative account of nonanalyticities in the finite-system en- 
tropy as observed for the special cases in [2, |3, |4[ ■ Then, 
inspired by a calculation in the density of states is 
split into two terms: The first is a sum of the leading 
nonanalytic contributions stemming from all the saddle 
points, whereas the second contains the "harmless" rest. 
Performing the thermodynamic limit of the thus obtained 
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expression, one can understand how, and under which 
conditions, the distribution of saddle points (as a func- 
tion of the potential energy per particle) may cause a 
phase transition to take place. These results can be inter- 
preted in terms of a recent approach relating phase tran- 
sitions and configuration space topology, and we apply 
our findings to the spherical model, explaining the puz- 
zling behavior of certain topological quantities reported 
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Preliminaries. — We consider classical systems of N 
degrees of freedom, characterized by a Hamiltonian func- 
tion of standard form, 



for any e > 0, with the universal function 



H(p,q) 



1 N 

It 



Pi + V(q), 



(1) 



where p — (pi,...,pjv) is the vector of momenta and 
q = (qi, . . . , qisr) the vector of position coordinates. The 
restriction to a quadratic form in the momenta is not 
essential, but simplifies the following discussion. The po- 
tential V is an analytic mapping from the configuration 
space IV C 'SL N onto the reals, and in general V will 
have a number of critical points (or saddle points) q c 
with vanishing exact differential, dV(q c ) = 0. We will 
assume in the following that all critical points of V have 
a non-singular Hessian Sjy- In this case, V is called a 
Morse function, and one can argue that the restriction 
of the potential to this class of functions is not a serious 
limitation: Morse functions form a dense, open subset of 
all analytic functions, and hence are generic. Therefore, 
even if V is not a Morse function, it can be made into 
one by adding an arbitrarily small perturbation 15]. 

Finite system nonanalyticities. — Our aim is now to 
investigate the analyticity properties of the (configura- 
tional) density of states 



(2) 



Q N (v) = / dqS(V{q) - Nv) 



or, equivalently, of the entropy sn(v) = ln[Ojv(u)]/iV of 
a classical system of N degrees of freedom, characterized 
by a Morse function V. The Morse property guarantees 
that all critical points of V are isolated. We thus may 
reduce the discussion to the effect of a single critical point 
of V on the analyticity properties of the density of states 
£In(v) as a function of the potential energy per degree of 
freedom v. 

Proposition 1 Let V : G — > R be a Morse function with 
a single critical point q c in an open region G C M. N . 
Without loss of generality, we assume V(q c ) = 0. Then 
there exists a polynomial P of degree less than N/2 such 
that at v = the density of states can be written in the 
form 



n N (y) = P(v) 



l>N.k 



(V) 



y/\det [S) V (q c 



+ o{v N ' 2 ^) (3) 



h N .k{v) 



(Ntt) n / 2 

r(f) 



{-l) — (-v)—Q(-v) forN,kodd, 

k N-2 



x < (-l)i«-T- G(v) 



for k even, 



(4) 



fc+l N—2 



(— 1) 2 v 2 7r 1 In for N even, k odd. 



Here k denotes the index of the critical point and is 
the Heaviside step function. 

This proposition gives a complete account of all types 
of nonanalyticities which can occur in fijv as the con- 
sequence of a single critical point of a Morse function 
V. In the presence of several critical points, their non- 
analytic contributions simply have to be added. For a 
proof of Proposition [l] the density of states O^r is calcu- 
lated separately below and above the critical value v = 0. 
By complex continuation it is possible to subtract both 
contributions and to evaluate the leading order of the 
difference 16 1. 

The index k of a critical point q c is the number 
of negative eigenvalues of the Hessian Sjv at q c , and 
it determines whether a critical point is a minimum 
(k = 0), a maximum (fc = N), or a proper saddle 
(k G {1, . . . ,N— 1}). Remarkably, the type of nonanalyt- 
icity described by hj^.k in (TJ|) does not depend on the pre- 
cise value of k, but only on whether iV and k are odd or 
even. One can verify that in any of the three cases in Q , 
fijv is I ^j^J -times continuously differentiable. This re- 
sult is in agreement with the nonanalytic behavior of the 
exact solution for the density of states of the mean-field 
spherical model as reported in 0. In other words, the 
density of states £In becomes "smoother" with increas- 
ing number of degrees of freedom, and already for moder- 
ate N it supposedly will be impossible to observe such a 
finite-system nonanalyticity from noisy experimental or 
numerical data. At first sight one might therefore suspect 
that the nonanalyticities of the entropy are irrelevant for 
large systems and have no effect in the thermodynamic 
limit, but the following considerations will show that such 
an assertion is premature. 

Thermodynamic limit. — In a recent letter [5j , Franzosi 
and Pettini have discussed for a certain class of short- 
range models the relation between nonanalyticities of the 
entropy or the free energy in the thermodynamic limit 
and topology changes of the subsets 



M v = {qeT N \ V(q) < Nv} 



(5) 



of configuration space Tn . Loosely speaking they found 
that, in order to have a phase transition at some poten- 
tial energy vt, a topology change of M v is necessary to 
occur at v — v t . Furthermore, in the Euler characteristic 
of M v of several models analyzed, a signature was found 
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at the phase transition energy t> t , corroborating the exis- 
tence of a relation between phase transitions and topol- 
ogy changes H, H| . It is a central proposition of Morse 
theory that, for potentials V having the Morse property, 
each topology change of M. v at some value v = v c corre- 
sponds to one or several critical points q c of V with crit- 
ical value v c — V(q c )/N [lfj. Combining these results 
with Proposition [U we infer that, despite their decreas- 
ing strength for large system sizes N, nonanalyticities of 
the finite-system entropy appear to be related to their 
infinite-system counterparts in some way. 

This relation, however, is not one-to-one: critical 
points of V are necessary, but by no means sufficient for 
a phase transition to occur. For several models studied, 
the number of critical levels v c was found to increase un- 
boundedly with the number of degrees of freedom of the 
system, and the levels become dense on some interval in 
the thermodynamic limit 0, HI, Q • It is a pertinent open 
question to understand how, and under which conditions, 
nonanalyticities of the finite-system entropy may give rise 
to a phase transition in the thermodynamic limit. This 
issue which we will address in the following, and — since a 
single critical point of V is too "weak" to produce a phase 
transition in the thermodynamic limit — it is the distri- 
bution function of the very many, densely lying critical 
levels which will play important roles in this discussion. 

By virtue of PropositionQ] we split the density of states 
of a system with N degrees of freedom into two terms, 



n N (v) = A N (v) +B N (v), 



(6) 



where Bn contains the leading nonanalytic contributions 
from the critical points of V, 



Bn{v) 



h N ,k{v - v c ) 



«o ?cK) vl det [fiv(q c j] 



(7) 



The sums extend over all critical values v c of V and all 
critical values q c (v c ) corresponding to the respective crit- 
ical value. Due to the alternating signs of h,N,k in Q we 
may consider Bn as a weighted alternating sum of the 
critical points. Then it is maybe not too surprising that 
the Euler characteristic of M. v , being an unweighted al- 
ternating sum of the critical points, was found to signal 
the presence of a phase transition for some models. 



In order to perform the thermodynamic limit, we have 
to consider the entropy 

s(v)= lim ^-ln[A N (v)+B N (v)], (8) 

since, at least for a system with short-range interactions, 
we may expect this quantity to exist [111 ] . It is important 
to note that Bm is of order e N for N — > oo and therefore 
gives a finite contribution to s. Now we would like to 
deduce the (non)analyticity of the entropy s from the 
properties of the thermodynamic limit expressions 



a(y) = lim —hiA N (v), b{v) = lim —bxB N (v). 

N—>oo IX N—>oo IX 



To this end it is instructive to rewrite Eq. © as 
s(v) = max {a(v), b(v)} . 



(9) 



(10) 



As a consequence of this expression, we cannot expect 
to draw any conclusions on s from the knowledge of one 
single function a or 6, but only from the interplay of 
both. If, for example, we find a nonanalyticity in 6, this 
nonanalytic behavior may be visible in s in one case, 
but it may simply be overruled by a larger a in another 
instance. Furthermore, a nonanalyticity in s may arise 
from a crossover between a and b when the dominant 
part in the maximum changes from a to b (or vice versa) . 

A "density splitting" in [6(, similar to our Eq. ([6]), in- 
dicates that one may expect the analytic part in An 
to converge uniformly to an analytic function for well- 
behaved short-range potentials, whereas possible non- 
analytic contributions are dominated by B^. Further- 
more it appears plausible from the discussion in the pre- 
ceding paragraphs that, at least for short-range models 
as covered by the theorem in [5( , nonanalyticities of the 
function &, stemming from the nonanalytic contributions 
of the critical points of the potential V, are crucial for the 
occurrence of a phase transition in the thermodynamic 
limit. 

When performing the limit JV-^oowe can for simplic- 
ity restrict the discussion to the subsequence with N odd 
(N even can be treated in a similar way). The sum over 
the critical values v c in ([7]) is converted into an integral, 
yielding 



B N (v) (27re) 



N/2 



dv c Wo(v c ) - Af 2 (v c )] (V - V c f N -^' 2 + I d« c [Af 3 (v c ) - M(«c)] (»c - «) (Ar - 2)/2 

I 



(111 



where the Afj are the limit distribution functions of 
the critical points with index k — j (mod 4) where, 



as a consequence of (JT]), each point is weighted by 
IdetffjvKoWjl ■ Note that B^ depends on the differ- 



4 



ences between the densities of certain types of critical 
points, and to illustrate the implications of this fact, we 
will now apply Eq. (TlTj) to the spherical model. This will 
allow us to understand the puzzling observations made 
by Risau-Gusman et at 01 when studying the topology 
of the subsets M v of this model. 

Application to the spherical model. — This model was 
introduced by Berlin and Kac [l2| as an exactly solvable 
caricature of the Ising model of a ferromagnet. Its con- 
figuration space Tn is an (N — l)-sphere with radius y/N, 
and the potential reads 



JV 



V : r 



JV 



JijQiQj- 



(12) 



The degrees of freedom qi of the model are placed on a 
d-dimensional cubic lattice, and we consider a coupling 
matrix J with elements Jy = \ when i and j are neigh- 
boring sites on the lattice, and Jy = otherwise. The 
spherical model is solvable in the thermodynamic limit 
for arbitrary d, and a phase transition from a ferromag- 
netic phase at low temperatures to a paramagnetic phase 
at high temperatures (or energies) occurs for all d ^ 3. 

In a recent letter [7], Risau-Gusman et al. reported 
a study of the topology of the subsets M v Q of 
this model. Computing the deformation retract of the 
Mv, they found that their topology is determined by the 
spectral density of the coupling matrix J. In the limit 
N — > oo, this spectral density is given by 



c(v) = - 

7T 



dx cos(xv)[Jq(x)] c 



(13) 



where Jo is the Bessel function of order zero. Depend- 
ing on the dimension d, the function c has one or sev- 
eral nonanalytic points, however — inconsistent with an 
expected relation between topology changes and phase 
transitions — none of these points coincides with the tran- 
sition energy vt of the spherical model. 

Instead of following the approach in we analyze 
the subsets A4 V by means of Morse theory. Although 
the potential V of the spherical model in ([12")) is not a 
Morse function, it can be made into one by adding a 
small coupling to a magnetic field h. In the simultaneous 
limit h — > and N — > oo one can relate the critical 
densities TV} in (|11|) to the limit distribution c, yielding 
TVo = M = W 2 = W 3 . We thus find b = for the 
spherical model, which implies that the contribution of 
the critical points of V to the entropy vanishes in the 
thermodynamic limit. This explains why no effect of the 
nonanalytic points of c is visible in the thermodynamic 
functions of the spherical model. 

At first sight this finding seems in conflict with the 
previous discussion where the relevance of the 6-term was 
emphasized. This apparent contradiction can be under- 
stood by noting that the spherical model is not genuinely 



short-range. Although the interaction term in (jT2j) is re- 
stricted to pairs of nearest neighbors, the spherical con- 
figuration space renders the interaction effectively long- 
range. For long-range systems however, as discussed in 
[l3j |. the connection between phase transitions and con- 
figuration space topology is not valid in general and we 
cannot assume the function a in (jlOD to be analytic. Sim- 
ilar to the observations made in [13j for a different long- 
range model, we suspect the entropy s(v,m) as a function 
of the potential energy v and the magnetization to to be 
a nonconcave function, such that the nonanalyticity in 
s(v) = max m s(w,m) arises from the maximization over 
rn. 

Summary. — We have analyzed the relation between 
saddle points of the potential energy V of classical N- 
particle systems and the analyticity properties of ther- 
modynamic functions. For finite systems, each saddle 
point q c was found to cause a nonanalyticity in the en- 
tropy sn(v) at the value v — V(q c )/N of the potential 
energy, and the functional form of the nonanalytic term 
is specified in Proposition [TJ Since the number of sad- 
dle points is expected to grow exponentially with N for 
generic potentials, we arrive at the remarkable conclusion 
that the finite-system entropy may have a large number 
of nonanalytic points. Considering the contribution of 
very many saddle points becoming dense in the thermo- 
dynamic limit, we discussed how, despite the increasing 
differentiability of s^r, a continuous distribution of sad- 
dle points may lead to a nonanalyticity in the infinite- 
system entropy. Interpreting our findings in the spirit of 
the topological approach to phase transition, our results 
indicate under which conditions topology changes of the 
subsets M v may lead to a phase transition in the ther- 
modynamic limit. An application of these findings to the 
spherical model allowed us to understand the puzzling 
observations of Risau-Gusman et al. that, for this model, 
topological signatures and phase transition energy do not 
coincide. 

Conceptually, our study explores the connection be- 
tween phase transition theory and the study of energy 
landscapes, a rapidly developing field with applications, 
among others, to clusters, biomolecules, and glassy sys- 
tems 14| . One might hope to profit from the considerable 



knowledge on the relation between energy landscapes and 
dynamical properties in future work. 
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